Probing possible decoherence effects in atmospheric neutrino oscillations 



O 

o 
o 

(N 
£ 

Q\ 
> 

m 
in 
o 

(N 
O 

o 
o 

i 

(D 
43 



X 



E. Lisi a , A. Marrone a , and D. Montanino fc 
a Dipartimento di Fisica and Sezione INFN di Bari, Via Amendola 173, 1-70126 Bari, Italy 
b Dipartimento di Scienza dei Materiali dell'Universitd di Lecce, Via Arnesano, 1-73100 Lecce, Italy 

(February 1, 2008) 

It is shown that the results of the Super-Kamiokande atmospheric neutrino experiment, interpreted 
in terms of <-> v T flavor transitions, can probe possible decoherence effects induced by new physics 
(e.g., by quantum gravity) with high sensitivity, supplementing current laboratory tests based on 
kaon oscillations and on neutron interferometry. By varying the (unknown) energy dependence 
of such effects, one can either obtain strong limits on their amplitude, or use them to find an 
unconventional solution to the atmospheric v anomaly based solely on decoherence. 



PACS: 14.60.Pq, 04.60.-m 

The Super-Kamiokande (SK) atmospheric neutrino 
experiment has found convincing evidence [jj] for the 
quantum-mechanical phenomenon of v flavor oscillations 
Q in the <-> v T channel. Such evidence consistently 
emerges from different SK data samples (sub-GeV lep- 
tons, multi-GeV leptons, and upgoing muons ||) as well 
as from other atmospheric v experiments Q . 

The simplest model for <-> v r oscillations involves 
two neutrino states v\ = (1,0) T and — (0, 1) T with 
masses m\ and m%, and two flavor states = (cg,s$) T 
and v r — (—sg,ce) T , where 9 is the neutrino mixing an- 
gle, c = cos, s = sin, and T denotes the transpose. The 
Liouville equation for the v density matrix p, 



P : 



-i\H, 



(1) 



is then governed (in the mass basis) by the Hamiltonian 
H = | diag(-fc, +k), where k = Am 2 /2E, Am 2 = ml - 
m\, and E (3> mi^) is the v energy (in natural units). 
The solution pit) of Eq. (Q) , with initial conditions p(0) = 
II„ fi (where H V/1 = v^®v^ is the f M state projector), gives 
the survival probability after a length x(~ t), 



is 2 n 



coskx) , (2) 



which is the well-known oscillation formula . 

Equation (^) beautifully fits the SK data H over a 
wide range of v energies (E ~ 10 _1 -10 3 GeV) and flight 
lengths (x ~ lOMO 4 km), provided that Am 2 ~ 3x 1CT 3 
eV 2 and s\ e ~ 1 |,§. Such striking agreement severely 
constrains possible deviations from the standard hamil- 
tonian H 1^,0. In this work we show that the SK data 
can also be used to probe deviations from the standard 
Liouville dynamics in Eq. ([!]), that might be induced by 
new physics beyond the standard electroweak model. 

In general, modifications of Eq. ([[]) emerge from dissi- 
pative interactions with an environment @], and can be 
parametrized by introducing an extra term D[p], 



p=-i[H,p]-V[p] , 



(3) 



which violates the conservation of Tr(p 2 ) and allows tran- 
sitions from pure to mixed states. The operator V has 



the dimension of an energy, and its inverse defines the 
typical (coherence) length after which the system gets 
mixed jjj. 

Among the possible sources of decoherence, a particu- 
larly intriguing one might be provided by quantum grav- 
ity, as suggested by Hawking in the context of black- hole 
thermodynamics fl(i|| . From such a viewpoint, any phys- 
ical system is inherently "open," due to its unavoidable, 
decohering interactions with a pervasive "environment" 
(the spacetime and its Planck-scale dynamics Fol- 
lowing the pioneering paper |L2"| , quantum gravity deco- 
herence effects have been investigated in oscillating sys- 
tems which propagate over macroscopic distances (see 
]i~3| for reviews). Analyses have been mainly focused on 
KK oscillations |l^,[l4],[L5| and on neutron interferometry 
p2] , ^6| , by assuming reasonable phenomenological forms 
for T>. In both systems, no evidence has been found for 
T> ^ 0, and strong limits have been derived on the quan- 
tities parametrizing V Jl3|: 



V\ \ < 1(T 21 GeV (KK, n systems) 



(4) 



Theoretical estimates for are very uncertain [Q\ 2[ , 
and can range from unobservably small values up to 
the limits in (Q). Therefore, it is wise to adopt a phe- 
nomenological viewpoint, trying to learn from experi- 
ments and to improve the laboratory limits (|j) with novel 
approaches, such as those provided by v oscillations. In- 
deed, attempts have been made to explain the solar v 
deficit through decoherence [pr|dT9| . It has also been sug- 
gested that decoherence might play a role in interpreting 
the atmospheric v data (]1]]2(| although, to our knowl- 
edge, no detailed analysis of the SK results has been 
attempted so far. The crucial point is that, for typi- 
cal atmospheric v energies (10 0±1 GeV), the oscillation 
length A = 2-njk spans the range ~ 10 3±1 km; then, if 
the (de)coherence length I is of comparable size, terms 
as small as ||-D|| ~ £~ x ~ 10~ 22±1 GeV can be probed. 

In order to fix a well-defined framework, we specialize 
Eq. (||) under reasonable (although not compelling) phe- 
nomenological assumptions. The most general require- 
ment is perhaps that of complete positivity f2UE2], cor- 
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responding to assume a linear, Markovian, and trace- 
preserving map p(0) — » This implies the so-called 
Lindblad form |23| for the decoherence term, 

V[p] = EJp, D n Dl) - 2D nP Dt , (5) 

where the operators D n arise from tracing away the en- 
vironment dynamics (see p4j for a recent proof). Master 
equations of the Lindblad form are ubiquitous in physics 
(see ||^5| for theorems and applications). Concerning 
v oscillations, such equations describe v propagation in 
dissipative media as, e. g., m atter with fluctuating density 
p6| or thermal baths P7| . Here, however, the environ- 
ment embeds possible new physics (e.g., the spacetime 
"foam" |ll[]) for which there is no established theory. 

In the absence of first-principle calculations, we assume 
that at least the laws of thermodynamics hold in the v 
system. The time increase of the von Neumann entropy 
S(p) = — Tr(plnp) can be enforced by taking D n = D\ 
||, so that Eq. (|) becomes T>[p] = £„[£>„, [D n , p}]. 
The conservation of the average value of the energy 
[Ti(pH)] requires, in addition, that [H,D n ] = @,[||. 

The hermitian operators p, H v , H, and D n , can be 
expanded |^] onto the basis formed by the unit matrix 
1 and by the Pauli matrix vector cr = (eri, cr 2 , 03 ) T . We 
take p = |(1 + p ■ cr), IL, M = i(l + q ■ cr), H = ~fc ■ cr, 
and D n = \d n ■ cr, where q = (s2e,0, C2e) T and k = 
(0, 0, -k) T . Defining G = £„ K| 2 1 - d n ® d T n , Eq. (§) 
is transformed into a Bloch equation, p = fcxp-G-p, 
which has a simple physical interpretation: the standard 
term k x p induces v oscillations, while the decoherence 
term G ■ p is responsible for their damping |^J2^] . 

The requirement [H, D n ] = implies that each vector 
d n is parallel to k |2^]. Therefore, the tensor G takes the 
form G = diag( 7 , 7 ,0) with 7 = E„ Kl 2 > & Thc 
general solution [p(<) = V" • p(0)] of the Bloch equation 
is then given by the evolution operator 

/ +er'< t cos kt +e" 7t sinfci 0\ 
V= -e-^sinkt +e-T*costo . (6) 

V 01/ 

If the system is prepared in the pure (zero entropy) 
state \p(0) = q], the asymptotic final state is p{oo) — 
(0,0,020). Since Tr[p 2 (oo)] = (1 + c% g )/2 < 1 and 
S^oo)] = — Cglncg — Sglnsg > 0, the system evolves 
indeed into a mixed state with positive entropy. Maxi- 
mal entropy (S = In 2) corresponds to maximal v mixing 
(s\q = 1). Purity and entropy are conserved only if p is 
prepared in a pure mass eigenstate [p(0) = (0, 0, ±1) T ]. 

The survival probability P M/J = i(l + q T ■ V ■ q) reads 

P w = 1 - §4,(1 - e-^ cos kx) , (7) 

which reduces to the standard expression (|[) in the limit 
7^0. For 7a; ~ O(l), one expects significant deviations 
from the usual oscillation fit to the SK data. 



We make a quantitative study of the effects of P MM 
in (0), by computing the theoretical SK lepton distribu- 
tions in zenith angle (d), and by fitting them to the SK 
data through a \ 2 statistics, as extensively discussed in 
||. The main difference from || is: (i) the 30 data bins 
for the SK distributions refer to a longer detector expo- 
sure (52 kton-year ||]); (ii) the oscillation probability is 
here taken from Eq. ([?]). In thc fit, we study both the 
case with (Am 2 , s\ 9 ,~f) unconstrained (oscillations plus 
decoherence) and the case with Am 2 = and (s 2 e , 7) un- 
constrained (decoherence only). We find significant dif- 
ferences in the results, depending on the energy variation 
assumed for 7 (which is not necessarily a constant param- 
eter). For definiteness, we discuss only three scenarios, 
corresponding to a possible power-law dependence of the 
kind 7 = 7o(£/GeV)™ with n = 0, 2, and -1. 

For n — (7 = 70 = const) the best fit with oscillations 
plus decoherence (x 2 nin = 22.6) is reached for Am 2 = 
3 x 10~ 3 eV 2 , s\ e — 1, and 70 = 0, which corresponds to 
the case of pure <-* v T oscillations. Since no evidence is 
seen to emerge for decoherence effects, meaningful upper 
bounds on the parameter 7 can be placed. By taking 
X 2 — Xmin = 6.25 (corresponding to 90% C.L. for three 
degrees of freedom) , we get 

70 < 3.5 x 10~ 23 GeV (n = 0) . (8) 

The limits at 95% and 99% C.L. are found to be 4.1 x 
10~ 23 GeV and 5.5 x 10" 23 GeV, respectively. The bound 
(||) shows that: (i) if decoherence effects have the same 
origin (e.g., quantum gravity) and similar size in the dif- 
ferent K, n, and v systems, then atmospheric v observa- 
tions can improve the current laboratory limits (^); and 
(ii) decoherence effects, if any, can develop only over a 
typical length scale £ — > 5600 km. 

Figure [l] shows (for n = 0) the zenith distributions of 
SK events for best-fit standard oscillations (70 = 0) and 
in the presence of an additional decoherence term (70 = 
10~ 22 GeV). The electron (e) distributions are unaffected 
(P ee = 1). In the sub-GeV p sample, decoherence is 
almost unobservable, due to the large intrinsic smearing 
§ of both energy and angle. In the multi-GeV p sample, 
the transition from no oscillation (P w ~ 1 for cost? ~ 1) 
to averaged oscillations (P MM ~ 1/2 for cost? ~ —1) is 
made only slightly faster by decoherence effects. Such 
effects are instead dominant in the higher-energy sample 
of upgoing p, where the oscillation phase kx is small, 
and decoherence generates a much faster suppression of 
vertical muons (cost? ~ — 1), corresponding to the longest 
v flight lengths. Finally, we find a bad fit (x 2 > 49) 
when oscillations are switched off [H = 0, corresponding 
to k = in Eq. (jTj)]. Therefore, in the case n = 0, the 
SK data cannot be explained solely by decoherence. 

The case n — 2 may also be of phenomenological in- 
terest, in the light of a possible dimensional guess of the 
form 7 cx E 2 /Mp pl| . In this case, decoherence effects 
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are even more disfavored than for n — 0, since they pro- 
duce a faster suppression of muons with increasing en- 
ergy, contrary to observations. We find an upper limit 
70 < 0.9 x 1CT 27 GeV at 90% C.L. [to be compared with 
the limit (@)]. For k — (decoherence without oscilla- 
tions) the fit is also very bad (x 2 > 70). 



Super-Kamiokande (52 kTy) 
zenith distributions of 
v events for 7=-y =const 
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FIG. 1. Effects of decoherence (70 7^ 0) on the distributions 
of lepton events as a function of the zenith angle (■&). The SK 
data are shown as dots with ±lcr error bars. The histograms 
represent our theoretical calculations. In each bin, the elec- 
tron (e) and muon (fi) rates R are normalized to standard (no 
oscillation, no decoherence) expectations Rq. 

From the previous cases (n = and n — 2) we learn 
that decoherence effects can be strongly constrained, the 
more the faster they increase with energy. Conversely, 
we expect weaker constraints for a decreasing energy de- 
pendence, such as for 7 oc E~ x (n = —I). 

The case n = — 1 may also be motivated by assum- 
ing that the exponent in Eq. ([?]) behaves as a Lorentz 
scalar. A boost from the v rest frame to the laboratory 
frame would then introduce a factor m v /E (just as for 
the oscillation phase), giving a decoherence parameter 
of the form 7 = 7o(-E/GeV) _1 . Of course, this ansatz 
should be taken with a grain of salt, since dissipative 
equations are known to entail problems with Lorentz in- 
variance 14 3^] (however, see p9| , |33| |). In any case, as- 
suming 7 = 7o(£'/GeV) _1 , we have performed a fit to the 
SK data with (A 

m ^ 1 s 26 ' To) unconstrained. The best fit 
is reached, once again, for 70 = 0, but the upper limit on 
70 is now relatively weak, 70 < 2 x 10~ 21 GeV at 90% C.L. 
Therefore, for n = — 1, one may add sizable decoherence 
effects to oscillations, without destroying the agreement 
with SK data. Can one switch off completely oscilla- 
tions, and explain the data as a pure decoherence effect? 
The answer is surprisingly positive. For Am 2 — 0, the 
best agreement with the data is reached at s 2 ,^ = 1 and 
70 = 1.2 x 10- 21 GeV, with Xmi„/^DF = 27.1/(30 - 2), 
giving a good fit. This case represents a novel solution to 
the atmospheric v anomaly, based solely on decoherence. 

Figure ||show such "exotic" best fit (decoherence with- 
out oscillations) as compared to the "canonical" best fit 
(oscillations without decoherence). The two cases appear 



to be almost indistinguishable within errors, although 
they entail completely different physics. It is amusing to 
notice that, for the two best-fit cases of Fig. |[ the 
survival probability approximately read 



P„ 



|[1 
|[1 



■cos(+P -L/E)] 
exp(-/3 ■ L/E)] 



(pure oscillations), (9) 
(pure decoherence), (10) 



where E is in GeV, L is the v pathlength (km), and 
(3 ~ 7 x 10~ 3 GcV/km. Both cases have the same asymp- 
totic behavior, namely, (P^) ~ 1 (i) for small (large) 
L/E. For intermediate values of L/E, the strong differ- 
ences between the oscillating cosine factor and the mono- 
tonic exponential damping appear to be effectively sup- 
pressed by the large smearing in the v energy and an- 
gle, due to the interaction and detection processes in SK. 
Therefore, future long-baseline accelerator experiments 
(such as K2K, MINOS, and the CERN-to-Gran Sasso 
project 1 34 1) will be crucial in discriminating the above 



two functional forms for P^ , by revealing the oscillation 
(or damping) pattern now hidden by smearing effects. 
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FIG. 2. Comparison of best-fit scenarios for pure oscilla- 
tions (solid line, as in Fig. |l|) and for pure decoherence with 
7 oc 1/E (dashed line). 

Finally, we test the best-fit decoherence case of Fig. || 
against the negative results of current i/n — > v T appear- 
ance searches @. In the CHORUS and NOMAD ex- 
periments JIH one has {L/E) ~ 0.025 km/GeV and 
P flT = l- P w ~ \P{L/E) (for Am 2 = and s 2 ^ 1). 
Then the experimental limit P^ T < 1.3 x 10~ 4 |34j im- 
plies the upper bound f3 < 1.1 x 10 -2 GeV/km, which is 
compatible with the best-fit value j3 ~ 7 x 10~ 3 GeV/km. 

In conclusion, we have performed a phenomenological 
analysis of modifications of the Liouville dynamics, in 
the context of atmospheric <-» v T transitions. Within 
a simple model embedding the relevant physics (oscilla- 
tions plus decoherence), we have found that the Super- 
Kamiokande data can be a sensitive probe of decoher- 
ence effects (e.g., originated by quantum gravity), sup- 
plementing current laboratory tests based on K and n in- 
terferometry. Depending on the energy behavior assumed 
for such effects, one can either constrain them strongly, or 
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use them to explain the atmospheric v anomaly without 
oscillations. 

We thank J. Ellis, R. Floreanini, S. Pakvasa, and 
S. Pascazio for useful comments. 

Note added. After submission of this Letter, two 
related works appeared |}(|. We also noted a recent 
preprint |3?J suggesting an exceedingly small theoreti- 
cal estimate for 7 (~ k 2 /Mp), which would discour- 
age current experimental tests with neutrinos (as well 
as with kaons and neutrons). It seems to us that such 
estimate p7[ , being essentially based on a dimensional 
guess, should be presently considered with great caution. 
In the absence of both a full dynamical theory and of 
ab initio calculations for decoherence effects, any current 
ansatz may prove to be wrong. This fact warrants phe- 
nomenological analyses as ours, whose results, inferred 
from experimental data, remain valid independently of 
(uncertain) guesses about the origin and the size of 7. 
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